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Abstract. - The Fractional Langevin Equation (FLE) describes a non-Markovian Generalized 
Brownian Motion with long time persistence (superdiffusion), or anti-persistence (subdiffusion) of 
both velocity-velocity correlations, and position increments. It presents a case of the Generalized 
Langevin Equation (GLE) with a singular power law memory kernel. We propose and numerically 
realize a numerically efficient and reliable Markovian embedding of this superdiffusive GLE, which 
accurately approximates the FLE over many, about r = iVlgfe — 2, time decades, where N denotes 
the number of exponentials used to approximate the power law kernel, and b > 1 is a scaling 
parameter for the hierarchy of relaxation constants leading to this power law. Besides its relation 
to the FLE, our approach presents an independent and very flexible route to model anomalous 
diffusion. Studying such a superdiffusion in tilted washboard potentials, we demonstrate the 
phenomenon of transient hyperdiffusion which emerges due to transient kinetic heating effects. 



' Introduction. — Anomalous diffusion and transport 
O processes possess a rich variety of applications spanning 

i ^ m any different research fields from plasma physics and 
nonlinear dynamical systems to condensed matter physics, 

O^l biophysics, epidemiology, and even quantitative finance 
JH3] ■ There are several very different theoretical ap- 
proaches to describe anomalous diffusion, from continuous 

QQtime random walks (CTRW) 0] including Levy flights and 

{Nl Levy walks [TJ[2] to the Generalized Langevin Equation 
(GLE) 5 20,. CTRW-based anomalous diffusion involves 

y—i Such unusual concepts as fractal time and subordination 

^^to a random clock which does not possess a finite mean 
period. In the continuous space limit it is often associ- 
£>ated with Fractional Fokker-Planck Equations (FFPEs) 

*k> EQ • The corresponding Langevin equations are local in 
vandom time [3T] and describe a doubly random process 
C^with infinite memory. These latter Langevin equations 
should not be confused [14] with the Fractional Langevin 
Equations (FLEs) [5I [TlTfHIIT7Il2"U] . Random time clocks 
without mean period entail the remarkable phenomenon 
of weak ergodicity breaking ■ The position, or ve- 

locity increments in this approach are independent in all 
basic models. Ergodicity can also be broken in other ap- 
proaches, e.g. in nonlinear Brownian motion |24j . Fur- 
thermore, GLEs including FLEs present a quite differ- 
ent approach incorporating long-time correlations or anti- 
correlations of the position increments, similar to the frac- 



tional Brownian motion (fBm) devised by Mandelbrot and 
van Ness [25] which can also be derived from a FLE [T4] . 
This is a benchmark feature [ToTPZO] of the latter approach. 
Another one is that this approach does not rely on the 
concept of random time with divergent mean clock period 
[TS] and is almost always ergodic, except for ballistic GLE 
diffusion HOlUM!]. 

Generally such non-Markovian processes are not char- 
acterized completely by a master equation for conditional 
probabilities, or Fokker-Planck equations |26j . However, 
the corresponding Langevin equations specify the stochas- 
tic process completely, containing all the information on 
trajectories. Phenomenologically, the GLE 
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formally presents a Newtonian equation of motion for a 
particle with mass m and coordinate x (we consider a one- 
dimensional model), subjected apart from a regular force, 
f(x,t) = —dV(x,t)/dx, to a zero-mean stochastic force 
(which adds energy to the Brownian particle) and 
a non-local in time frictional, or dissipative force (which 
takes off energy from the Brownian particle). Both pro- 
cesses are balanced at thermal equilibrium, i.e. Brownian 
motion never ceases and obeys the fluctuation-dissipation 
theorem, FDT. In ([T]), the dissipative force is assumed to 
have the form of a linear velocity-dependent friction with 
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its memory characterized by the integral kernel r){t). The 
FDT is obeyed, when the noise is Gaussian [57] and the 
memory kernel and the noise autocorrelation function are 
related by the fluctuation-dissipation relation (FDR) [5] 
reading 

(C(t)C(t'))=k B Tr 1 (\t-t'\) (2) 

with T being the environmental temperature. 

Apart from this phcnomcnological justification, the 
GLE can be derived microscopically from a Hamiltonian 
model involving coupling of the diffusing particle to a ther- 
mal bath of harmonic oscillators obeying initially canon- 
ical distribution at temperature T [H1IZ1 E] - Anomalous 
diffusion can be related to a power-law memory kernel 

= lsin(7ra/2)| r + 

7T/Z 

which we write with a short-time cutoff 1 /lu c correspond- 
ing to the largest frequency of the bath oscillators u> c , 
a > 0, and T(x) is the gamma- function. By eq. ©, 
for < a < 2 such r\ (t) corresponds in the singular limit 
lj c — > oo to a fractional Gaussian noise with the Hurst 
exponent H = 1 — a/2 [25]. Then the solution of the GLE 
yields a fractional Brownian motion in the limit m — > 
[14] . In the case of free diffusion [f(x,t) = 0] and for 
< a < 2 independently of u c the noise-averaged position 
variance a 2 (t) — (Ax 2 (t)) = ([x(t) — (x(t))] 2 ) grows with 
time asymptotically as <r 2 (t) ~ 2fcsTi Q /[r/ a T(l + a)] [§]. 
This corresponds to sub-diffusion in the case < a < 1 
(sub- linear growth, sub-Ohmic thermal bath), normal dif- 
fusion for a = 1 (linear growth, Ohmic bath), and su- 
perdiffusion for 1 < a < 2 (super-linear growth, super- 
Ohmic bath). In terms of the integral frictional strength, 
77(0) = J Q r){t')dt' , these behaviors are intuitively clear: 
subdiffusion corresponds to 77(0) — > 00, normal diffusion 
to 77(0) = const and supcrdiffusion to f)(0) — > 0. For 
a > 2, the free diffusion is always ballistic, cr 2 (t) ~ t 2 , 
and noncrgodic because the velocity autocorrelation func- 
tion (VACF) does not decay to zero. Therefore, the emer- 
gence of hypcrdiffusion cr 2 (t) ~ t x , with A > 2 within GLE 
model is rather surprising. All these results can be eas- 
ily obtained from a general expression for the stationary 
VACF, found for an arbitrary memory kernel first by Kubo 
[5] and reading K v (s) = IcbT '/[ms + fj(s)] in the Laplace 
space, by taking into account that the position variance 
is the twice-integrated VACF. One has to remark at this 
point, that the Laplace-transformed rj(t) is i)(s) = ry Q s Q_1 
in this model in the limit uj c — > 00. However, the inverse 
Laplace transform for t)(s) = r\ a s a ~ x only exists for a. < 1, 
i.e. for normal and subdiffusive friction. For 1 < a < 2 
and a finite w c , the kernel n (t) in eq. ([3]) starts from a 
positive part and then becomes negative so that its to- 
tal integral is zero, independently of uj c . When the cutoff 
frequency tends to infinity, the memory kernel becomes 
singular, starting from a positive singularity and being 
negative otherwise. The frictional term can be recast in 



this limit with the help of a fractional Riemann-Liouville 
derivative, see below. 

For example, a spherical particle of radius R moving 
with velocity v(t) = x(t) in an incompressible liquid of 
kinematic viscosity u and density p experiences a hydro- 
dynamic force [55] 

F v (t) = -2irpR 3 (^x+^x 

R V t V-oo Vt^F J 

This classical result due to Boussincsq and Basset [11] 
generalizes the well-known by Stokes. The first term 
yields a mass renormalization of the Brownian particle 
?n — y 771 -\- Am with Am = 2irpR 3 /3, which is assumed 
to be implicitly done. It is present also in the absence 
of dissipation, i.e. for fx — > 0. The second term cor- 
responds to the Stokes friction, and the third term is 
due to a finite relaxation time r r = R 2 j \i of the dis- 
turbed velocity field of the liquid. It reflects hydrody- 
namic memory. An interesting mathematical interpreta- 
tion of this term can be given within the formalism of frac- 
tional derivatives [BUI]. Namely, using the definition of 

the Riemann-Liouville fractional derivative, t Q D t f(t) := 
rJT^yjlt It - t'V, < 7 < 1, acting on some 

test function f(t), it can be recast in the form [TT] 

F ad {t) = ~7 la ^I)"" 1 x{t) (5) 

with 7] a = rji R / y/Jl and a = 3/2, where 771 = dtirRpp is 
the Stokes friction coefficient. The corresponding GLE 
was termed FLE 1 1 . Independently of this interpreta- 
tion it is known |29) to yield the famous power-law decay 
of the VACF, which was revealed in molecular dynamic 
simulations by Alder and Wainwright |30j . 

The presence of a normal Stokes friction term makes the 
corresponding diffusion asymptotically normal. However, 
anomalous superdiffusivc motion can emerge on a tran- 
sient time scale t < r r = R 2 / p, for light particles [TT]. If 
to neglect ad hoc the Stokes term and to set x(t) = for 
t < to = (i.e. the particle starts to move at t = 0) the 
corresponding superdiffusive FLE reads [51111]. 

-a-l .. . dV(x.t) ., . 

mx + 7) a D t x(t) + ^ ' ; = C(t) (6) 

with a = 1.5. This FLE serves as one of the basic models 
for supcrdiffusion with 1 < a < 2. It corresponds to a 
GLE with a singular memory kernel, which in mathemat- 
ical sense is a generalized function. As discussed above, 
this kernel starts from a positive singularity at t = and 
then is negative, decaying to zero in accordance with a 
power law t~ a , so that its total integral is zero. 

The presence of a nonlinear time-dependent force 
f(x,t) = —dV(x 7 t)/dx modifies this well-established pic- 
ture considerably. Since general analytical results are then 



p-2 



Markovian embedding of fractional supcrdiffusion 



scarce and most likely generally nonexistent, the reliabil- 
ity of numerical simulations is a key issue. In particu- 
lar, numerically tractable models can be obtained by ap- 
proximating the given power law memory kernel by a fi- 
nite sum of exponentials, see in [T3J[TBJ[T5]. By means 
of this approximation it is possible to represent the non- 
Markovian dynamics in the (x, v) plane as a projection of 
a fully Markovian dynamics in a hyperspace of dimension 
D = N + 2, where N is the number of exponentials in the 
approximation |T^jT6l[l8] . Then the central point is that 
one can propagate the corresponding Markovian dynam- 
ics locally in time by very reliable algorithms with a well 
controlled numerical precision and by increasing N one 
can approximate the FLE dynamics ever better. Surpris- 
ingly one finds that the practical embedding dimension D 
need not be large to achieve an excellent approximation 
within statistical errors of stochastic simulations. More- 
over, this approach can be used independently of the FLE 
with some advantages: (i) The stochastic propagation is 
local in time and can readily be continued beyond end 
point. There is no need to generate a realization of long- 
correlated noise £(t) for the whole time span of simula- 
tion fixed in advance (without a possibility to continue) 
and solving the intregro-differential equation numerically 
for any such noise realization. This dramatically saves 
computer memory and enables extremely long simulations 
with appreciably small statistical errors, (ii) The cor- 
responding multi-dimensional Fokker-Planck equation is 
known explicitly for arbitrary f(x,t) which can be used 
to develop an analytical theory for nonlinear dynamics. 

Model. — We generalize the modeling of viscoelastic 
subdiffusion in ref. [TB] to superdiffusive case, 1 < a < 2, 
and approximate the superdiffusive memory kernel as 
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The first singular term mimics the positive singularity in 
the memory kernel of the FLE . The sum of exponen- 
tials obeys a fractal scaling with negative weights and ap- 
proximates the power law decay [5] of this memory kernel; 
i.e., r)(t) oc —t~ a , for t > 0, so that 77(0) = 0. Choosing 
Vq = uj Ci the power law regime extends in this approxima- 
tion from a short time (high-frequency) cutoff, t; = oj~ , 
to a large time (small frequency) cutoff, Tu = nb N , where 
b is a scaling dilation parameter. Such a fit is known to ex- 
hibit logarithmic oscillations superimposed on the power 
law |2J. Their amplitude is, however, small and can be con- 
trolled by the choice of b. By adjusting b and N for a given 
a one can approximate t~ a over about r — Nlgb — 2 time 
decades between two time cutoffs, which are always physi- 
cally present, beyond FLE modeling. Fig. Q] illustrates the 
quality of the approximation of t with the parameters: 
u = 10 3 , b = 5, N = 13, and C a (b) = 1.78167. One can 
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Fig. 1: Approximation of the power law behav- 
ior of the friction kernel by a sum of exponentials 
Ca(*0£i=i (V&T exp (-iW^) eq. ©] with the 

parameters: N = 13, u Q = 10 3 , b = 5, and C a (b) = 1.78167. 



detect a good agreement over about r ~ 7 decades in time. 
Clearly, with decreasing b and increasing N, one can fur- 
ther improve and control the quality of the approximation 
[16] which should be consistent with statistical errors of 
Monte Carlo simulations to avoid unnecessary numerical 
load. 

Next, we introduce N auxiliary variables u% and the 
corresponding multi-dimensional Markovian dynamics in 
the hyperspace of dimension D = N + 2, 



N 



x{t) =v(t) 

mv(t) = - V'(x, - £ Ui(t) - Vov{t) + \Z2fc B T7jb6>(t) 

Ui{t) = - r)iv(t) - ViUi{t) + y/2k B TriiVi£i(t), (8) 

where v.i = v /b l , rji = C a {b)r] a v? /\T(l - a)\, for i 

1 Y . and i]o = ^iLiVi/^i- Furthermore, £o(t) and 

the £i(t) are N + 1 delta-correlated white Gaussian noise 
sources of zero-mean and unit intensity. N of them 
are totally uncorrclated, (£i(t)£j(t')) = o~ij8(t — t'), for 
i,j = 1,...,N, and for i = 0,j = 0. However, the noise 
£0 (t) is chosen as a weighted, normalized sum of the other 
independent noises, 



(9) 



In the limiting case N = 1, our present model yields 
the minimal 3-dimcnsional embedding of ballistic CLE su- 
pcrdiffusion developed in Refs. P3BI3II an d presents thus 
a generalization of this earlier model to the sub-ballistic 
case. Given the lower integral limit to = in the GLE (in- 
stead of minus infinity), the exact reduction requires that 
the Ui(0) are independently Gaussian distributed with 
zero-mean and variance (^(0)^(0)) = ksTijiSij to en- 
sure the stationarity of the noise ((t) and FDR (fJJ) for all 
times. 
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Furthermore, on the time scale t > Th, the diffusion be- 
comes ballistic in our modeling. However, Th becomes 
exponentially larger with increasing N, which ensures, 
that the practical embedding dimension can be reason- 
ably small (D = 15 in simulations below). When speak- 
ing about asymptotic behavior below we yet assume that 
t < Th, but choose Th so, that it cannot be reached numer- 
ically. This procedure is similar to the subdiffusive case 

Numerical Simulations. — Below, we investigate 
superdiffusion with a = 1.5 in a tilted washboard potential 
of the form V(x) = — Vq cos(2irx / xq)— Fx. For sake of con- 
venience we transform the equations above into dimension- 
less units by scaling time t in units of To = (r^/m) 1 ^ 01-2 ^, 
distance in xq, energy (which applies to Vq, ksT and 
Fxq) in AE = m(xo/ro) 2 and Ui in tuxo/t^. As a con- 
sequence of the time scaling, vq must be scaled in 1/ro, 
which implies, that r]o/m is scaled in 1/tq and rji/m in 
1/tq. We integrated the dimensionless equations with a 
standard stochastic Euler algorithm using a combination 
of Mersenne-Twister and Box-Muller algorithms to gen- 
erate the Gaussian random numbers. In each simulation 
an ensemble of 10 4 particles was propagated with a time 
step At — 10 -4 to achieve (weak) convergence of the en- 
semble averaged results. The end point of simulations is 
£max = 10 4 . We used the friction kernel parameters of 
the approximation shown in fig. Q] and distributed the 
particle velocities initially thermally with (v(0)) = and 
(v 2 (0)) = ksT jm = v\. All particles were set initially to 
the position a;(0) = 0. 

Free superdiffusion. We first test our method by com- 
parison of the numerical results for free superdiffusion, i.e., 
Vq = and F = 0, with the available analytical solution 
of FLE [H], 

(Ax 2 (t)) = 2v 2 t 2 E 2 ^ 3 [-(t/T ) 2 - a ] . (10) 

Here, E a ^(t) = J2^=o t n /T(an + /3) is the generalized 
Mittag-Lcffler function. Furthermore, for rj(t) in eq. (|7|) 

the Laplace-transformed variance cr 2 (s) is a rational func- 
tion which can also be inverted to the time domain. The 
agreement in fig. [5] among the analytical FLE result, the 
analytical result for the Markovian embedding and the 
simulation results is indeed very good. Initially, the dif- 
fusion is always ballistic, (Ax 2 (t)) ~ t 2 , turning over into 
the asymptotic behavior (Ax 2 (t)) ~ t 1 . 

Superdiffusion in a tilted periodic potential. Next, we 
study superdiffusion in a tilted periodic potential, for 
which no analytical solution is available. Similar to the 
free case, the superdiffusion starts again with ballistic dif- 
fusion and asymptotically again yields, (Ax 2 (t)) ~ t 1 - 5 , 
cf. fig. [3] for F = 10. However, on an intermedi- 
ate time scale a hyperdiffusive regime is developed, with 
(Ax 2 (t)} ~ t x , where A > 2. Such an puzzling intermedi- 
ate regime (note that free superdiffusion cannot be faster 
than ballistic within the GLE description) with a highly 
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Fig. 2: Comparison of the analytical FLE solution for the 
mean-squared displacement in eq. (|10[l with the analytical 
results for the Markovian embedding and the corresponding 
numerical results, a = 1.5. 



enhanced power-law dependence of the position variance 
was also found recently for ballistic superdiffusion, see |19j , 
and seems to be a generic feature of driven superdiffusion 
in non-confining potentials. Within the GLE description 
with a power law memory kernel and a quite different nu- 
merical approach (not based on Markovian embedding), 
such a hyperdiffusion was revealed first in ref. [T5]. How- 
ever, it was perceived there as an asymptotic regime prob- 
ably because of insufficiently long time propagation. It is 
clear from fig. [3]that for a sufficiently small bias F = 0.25 
we also cannot arrive at the asymptotic regime within two 
weeks of simulations. For larger bias strengths this is, 
however, possible. The mean velocity of diffusing par- 
ticles also grows with time and when the corresponding 
mean kinetic energy becomes so large that the periodic 
potential ceases to play a role, the diffusion attains asymp- 
totically the regime of free superdiffusion. However, on 
an intermediate time scale such a confined superdiffusion 
can be much faster than the free one. This surprising phe- 
nomenon somewhat resembles giant acceleration of normal 
diffusion in washboard potentials (35] , but has a distinctly 
different origin [19]. 

The intermediate behavior of the position mean-squared 
displacement (MSD) is a consequence of a transient heat- 
ing of the particles from the bath temperature T to the ki- 
netic temperature 7ki n , which we define |19j via the width 
of the velocity distribution ksT] i i n (t)/2 = m(Av 2 (t))/2, 
Av(t) = v(t) — (v(t)). Such a kinetic temperature should 
not be confused with the thermodynamic temperature. 
Nevertheless, it presents a very helpful concept character- 
izing the kinetic energy of the disordered motion, when the 
kinetic energy of the directed motion is subtracted from 
the mean kinetic energy, in accordance with the physical 
meaning of temperature. 

Fig. 2] shows the time evolution of the velocity 
variance for the diffusion in fig. [3] One can im- 
mediately see that the transient hyperdiffusion is re- 
lated to the transient kinetic heating. To explore 
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Fig. 3: Position mean square displacement of particles in a 
tilted washboard potential for several different values of the 
bias strength F and the fixed potential amplitude Vo = 1 and 
the bath temperature T = 1. The dotted line corresponds to 
the free case. 



this relation in more detail we notice that the po- 
sition MSD equals the twice integrated VACF, i.e., 

(Ax 2 (t)) = 2f*dt'f*'dt"{Av(t')Av(t")). Introducing a 
normalized VACF K v (t",t') = (Av(t')Av(t"))/(Av 2 (t")) 
one can rewrite the integrand as (Av(t')Av(t")) = 
(Av 2 (t"))K v (t" , t') and single out the evolution of the ve- 
locity variance. If K v (t' ,t' + r) would remain constant, 
or it would tend asymptotically to a constant with in- 
creasing r (indicating the breaking of crgodicity), then a 
power law increase of the kinetic temperature with time, 
Tki a (t) ~ t@ , would yield to a hypcrdiffusive law for the 
position MSD (Ax 2 (t)) ~ T kin {t)t 2 ~ t 2+f3 . This is what 
occurs for ballistic superdiffusion in tilted washboard po- 
tentials which is nonergodic [T5J. There, the kinetic tem- 
perature increases until it saturates at a large final value 
|19j . In contrast, the studied generalized Brownian mo- 
tion is ergodic. Its free motion VACF does decay to zero. 
It also decays in the tilted washboard potentials, see in 
fig- 03 for the case F = 1.5 in figs. 13141 However, this 
decay is enormously retarded in the transient regime, cf. 
fig- El This retardation explains qualitatively the quasi- 
noncrgodic origin of hypcrdiffusion also in the present 
case. Indeed, if K v (t' , i' + r) were a constant, then the hy- 
pcrdiffusion power law exponent would beA = 2 + /3«2.8 
for the velocity variance growth depicted in fig. 0] This is 
not much different from the observed value A ps 3 in fig. 
[3] However, very different from the case of ballistic mem- 
ory kernel [TS], the kinetic temperature, associated with 
the velocity variance, starts to decline towards the bath 
temperature, see in fig. 3J This is because the motion is 
ergodic and after the periodic potential ceases to play a 
role the velocity ACF has to decay to zero much faster (like 
in the free case). Since the dissipation is much stronger 
than in the ballistic case, the particles are not strictly ac- 
celerated by the bias F. Their mean velocity grows rather 
sublinearly, (v(t)) ~ Ft"^ 1 and one can show that for any 
a < 2 the mean kinetic energy of particles becomes neg- 
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Fig. 4: Mean square fluctuation of the particles velocity for the 
diffusion of fig. [3] 

ligible in the course of time as compared with the work 
done by the biasing force F on its increase. This is very 
different to the ballistic case, where asymptotically a fi- 
nite portion of work is used for the heating (see online 
Supplement of ref. [HI), and leads eventually to the de- 
cline of the kinetic energy of disordered motion, cooling 
back to the bath temperature T. If our explanation of 
the transient hypcrdiffusion in terms of a delayed decay 
of VACF is correct, it should be able to explain another 
feature in fig. [3] Namely, that the hypcrdiffusive regime 
first turns over into a transiently decelerated supcrdiffu- 
sion with A ps 1.2 < a (see for F = 1.5), before A grows 
back to a. Indeed, after reaching the maximum in fig. U 
the velocity variance starts to decay in accordance with a 
power law j3 ps —0.9. Then, our reasoning with a strongly 
delayed VACF decay yields A = 2 + /3ps1.1 which is close 
to the observed A ps 1.2 in fig. [3J This confirms that our 
line of reasoning is consistent. 

Summary and conclusions. — In this work we pro- 
posed a general and simple method for Markovian embed- 
ding of a supcrdiffusivc non-Markovian CLE dynamics and 
showed that it can be used to approximate FLE dynamics 
over many time decades, serving also as an independent 
approach to model supcrdiffusion. We studied numerically 
such a superdiffusion with a = 1.5 in a tilted washboard 
potential. Concordant with our prior findings for the bal- 
listic diffusion case [19] . a hypcrdiffusive transient regime 
was found. We gave a simple physical explanation of this 
transient regime as a kinetic heating effect in terms of the 
growing velocity variance and a strongly delayed decay of 
the velocity autocorrelation function. This is similar to the 
case of nonergodic ballistic supcrdiffusion in washboard 
potentials. The transient can be very long and supcrdiffu- 
sion can become enormously accelerated, compared to the 
free superdiffusion, during this transient regime. However, 
very different from the ballistic case, the transient heating 
is followed by a subsequent cooling back to the tempera- 
ture of the thermal bath, after the kinetic energy of the 
particles, which grows in time, exceeds much the potential 
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Fig. 5: Normalized velocity autocorrelation function for the 
free superdiffusion (dotted line, analytical result) and for the 
biased diffusion of fig. [3] for F = 1.5. One can see, that the 
latter decays much more slowly in the transient regime, for two 
values of to, which correspond to the start of the ascent of the 
velocity variance in fig. 21 to = 4, and the descent, to = 40. In 
the asymptotic regime, to = 9500, the VACF is again the same 
as in the free case. 

energy. This cooling effect is due to the fact, that the mo- 
tion remains ergodic in the studied case and the friction is 
sufficiently strong to take off the extra part of the kinetic 
energy, which was built up during the transient heating 
regime. During this transient cooling regime the power 
law exponent of the diffusion becomes less than the one of 
the free superdiffusion and gradually grows to the latter 
in the course of time. Then the VACF coincides asymp- 
totically with the one of the free motion. The influence of 
the periodic potential becomes forgotten, very differently 
from the nonergodic ballistic case. 

In conclusion, we expect that our general methodology 
will be used in a number of future applications of anoma- 
lous GLE diffusion. Moreover, the surprising transient 
heating/cooling effects are expected to attract a further 
attention not only of theorists but also of the experimen- 
tal community. 

* * * 
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